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ABSTRACT 



Context. X-ray surveys facilitate investigations of the environment of AGNs. Deep Chandra observations revealed that the AGNs source surface 
\^ ' density rises near clusters of galaxies. The natural extension of these works is the measurement of spatial clustering of AGNs around clusters 
and the investigation of relative biasing between active galactic nuclei and galaxies near clusters. 

Aims. The major aims of this work are to obtain a measurement of the correlation length of AGNs around clusters and a measure of the averaged 
clustering properties of a complete sample of AGNs in dense environments. 

Methods. We present the first measurement of the soft X-ray cluster-AGN cross-correlation function in redshift space using the data of the 
O ROSAT-l<iEP survey. The survey covers 9x9 deg^ around the North Ecliptic Pole where 442 X-ray sources were detected and almost completely 
^ spectroscopically identified. 

CZ3 Results. We detected a >3cr significant clustering signal on scales <50 h^^ Mpc. We performed a classical maximum-likelihood power-law 



Oh 



fit to the data and obtained a correlation length So=S.l*g^ h^^ Mpc and a slope 7=1.7;^[J-^ ( Icr errors). 

Conclusions. This is a strong evidence that AGNs are good tracers of the large scale structure of the Universe. Our data were compared to 
the results obtained by cross-correlating X-ray clusters and galaxies. We observe, with a large uncertainty, a similar behaviour of the AGNs 
clustering around clusters similar to the clustering of galaxies around clusters. 

Key words. Galaxies: clusters:general - Galaxies: active - X-rays: galaxies: clusters - Cosmology: large-scale structure of Universe - 
Cosmology: dark matter 



1. Introduction 

The current paradigm of galaxy formation assumes that all 
types of galaxies reside in dark matter (DM) haloes, and that 
the properties of these haloes determine to some extent the 
properties of the galaxies inside (White & Rees 1978). In gen- 
eral, the clustering amplitudes of haloes depend on halo mass. 
The relation between the clustering properties of both DM 
haloes and galaxies (biasing) should thus tell us something 
about the physical processes leading to the formation and evo- 
lution of galaxies. Clusters of galaxies are the highest peaks 
in the global mass distribution of the Universe and should fol- 
low a direct and simple biasing scheme - mainly related to the 
underlying primordial Gaussian random field (Kaiser 1987). A 
simple means towards a better understanding of galaxy bias- 
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ing is thus provided by studies of the relative biasing between 
galaxies and clusters of galaxies. 

As a first step, previous investigations estimated two-point 
statistics like the oMfo-correlation function (MuUis et al. 2004, 
GilH et al. 2005, Basilakos et al. 2005, Yang et al. 2006). They 
could show that, in fact, AGNs trace the underlying cosmic 
large-scale structure. In addition, the large-scale structure of 
X-ray selected galaxy clusters could be studied in some detail 
(e.g. Schuecker et al. 2001), but without investigating the link 
between clusters and AGNs. 

In the present paper, we concentrate on the study of the rel- 
ative clustering between X-ray selected AGNs and galaxy clus- 
ters. Our work improves on most previous work on the large- 
scale structure of X-ray selected AGNs in two important as- 
pects. First, with the exception of Mullis et al. (2004), our sam- 
ple is the only one that is spectroscopically complete (99.6%). 
GilH et al. (2005) used the CDFS ( 35%) and the CDFN (50%). 
The Basilakos et al. (2005) sample had almost no spectroscopic 
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redshifts. Yang et al. (2006) used the CLASXS sample (52% 
complete) and the CDFN (56% complete). 

Another motivation for our work is that over the last sev- 
eral years. X-ray observations revealed that a significant frac- 
tion of high-z clusters of galaxies show overdensities of AGNs 
in their outskirts (i.e. between 3 h^^ Mpc and 7 h^^ Mpc from 
the center of the cluster) (Henry et al., 1991; Cappi et al., 2001; 
Ruderman & Ebeling 2005, Cappelluti et al., 2005, and ref- 
erences therein). These overdensities were however detected in 
randomly selected archive targeted observations of galaxy clus- 
ters. While these overdensities are highly significant (up to 8cr) 
when compared to cluster-free fields, the incompleteness of the 
samples does not allow drawing any conclusion about the aver- 
age clustering properties of AGNs around clusters. The major- 
ity of the sources making these overdensities have no spectro- 
scopical identification and therefore any information on their 
spatial clustering is lost. More recently Branchesi et al. (2007) 
showed that at high-z the source surface density of AGNs sig- 
nificantly increases even in the central regions of the clusters. 
These results imply that further progress will come from study- 
ing the three dimensional spatial distribution of AGNs around 
clusters. A natural way to characterize this specific type of 
clustering is given by the three-dimensional croii-correlation 
of AGNs and galaxy clusters, the computation of which needs 
complete redshift information for all objects, which is rare in 
X-ray surveys. 

In this respect, the ROSAT North Ecliptic Pole (NEP) sur- 
vey (Henry et al., 2001,2006; Voges et al., 2001) is one of the 
few X-ray surveys covering a sufficiently large volume with an 
almost complete follow-up identification of AGNs and clusters 
(i.e. 440 sources spectroscopically identified of 442 detected). 
This survey thus provides a very useful basis for more precise 
investigations of the relative clustering properties of these two 
types of objects. 

We organized the present paper in the following way. In 
Sect. 2, we describe the ROSAT-NEP survey data which we 
use for our investigations of the spatial distribution of X-ray 
selected AGNs and galaxy clusters. For the statistical analy- 
sis we estimate their cross-correlation. A useful estimator for 
this statistic and the mock samples needed for its determi- 
nation are described in Sects. 3 and 4, respectively. The re- 
sults are presented in Sect. 5, and are discussed in Sect. 6. In 
this paper, we assume a (concordance) Friedmann-Lemaitre 
Universe characterized by the Hubble constant given in units 
of /170 = //o/(70kms"' Mpc"'), the normalized cosmic matter 
density = 0.3, and the normalized cosmological constant 
- 0.7. Unless otherwise stated, errors are reported at the 
Icr confidence level. 

2. The data 

The ROSAT NEP survey covers a region of 9 x 9 deg^ around 
the North Ecliptic Pole (17''15'" <a< 18''45'", 62° < 6 < 71°) 
observed with the PSPC proportional counter as part of the 
ROSAT All Sky Survey (Henry et al., 2001; Voges et al., 2001) 
with a flux limit of 2xl0"''*erg cm"^ s"' in the 0.5-2 keV 
energy band. 442 X-ray sources were detected and 440 opti- 
cally identified. Spectroscopic redshift information is available 



for 219 AGNs and 62 clusters of galaxies. The clusters have 
redshifts z <0.81 with a median of 0.18 and the AGNs have 
z <3.889 with a median of 0.4 (Fig.[T]l. For the purpose of this 
work we selected all the clusters and the 185 AGNs with z <1 
(Gioia et al. 2003). 

Such a dataset was used also by Mullis et al. (2004) for the cal- 
culation of the 3-D auto correlation function of X-ray selected 
AGNs. Mullis and collaborators find significant clustering on 
scales smaller than ~ 43 /i^,} Mpc with a correlation length of ~ 
10.4 /i^y Mpc, and a slope of the correlation best-fit power law 
of 7=1. 8. 

3. Cluster-AGN spatial cross-correlation 

The cross-correlation function fcA of clusters and AGNs is de- 
fined by the joint probability to find, at a distance r, one clus- 
ter in the infinitesimal comoving volume element 6Vc and one 
AGN in the comoving volume element 6Va, 



6P^ncnA[l+^CA(r)] 6Vc6Va 



(1) 



where «c and are the mean comoving number densities of 
clusters and AGNs, respectively. In calculating the differential 
cross-correlation in redshift space we used an adapted version 
of the Landy-Szalay estimator (Landy & Szalay, 1993; see also 
e.g. Blake et al., 2006), 



DcDa - RcDa - RaDc + RaRc 
RaRc 



(2) 



where DcDa, RcDa, RaDc and RaRc are the normalized num- 
ber of pairs in the /-th redshift space separation s bin for the 
clusters data-AGNs data, cluster random-AGNs data, AGNs 
random-clusters data and clusters random-AGNs random sam- 
ples, respectively. Using the symbols D and R to represent the 
data and random samples, respectively, and C and A to iden- 
tify clusters and AGNs, respectively, the normalized pairs are 
expressed by 



DcDa 



n 
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(Si) 



(Si) 
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npair.RcD,x(Si) 



Here, A^o, , Ndc . -^r,, , Nr^ are the total numbers of AGNs and 
clusters in the data and in the randomly generated samples, re- 
spectively. The quantities «pa„. represent, adopting the symbol- 
ism used above, the actual number of pairs measured in the 
random and data samples as a function of the redshift space 
separation i,. The distances were computed assuming for the 
position of the clusters the centroid of the X-ray emission while 
for AGNs the optical positions were adopted. In order to have 
a good signal-to-noise ratio (SNR) the data were grouped in 
logarithmic bins of Alog(i/!^g Mpc)=0.15. 



4. Random Samples 

Equation |2] indicates that an accurate estimate of the distribu- 
tion function of the random samples is crucial in order to ob- 
tain a reliable estimate of ^ca- Several effects must be taken 
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into account when generating a sample of objects in a flux lim- 
ited survey. Simulated AGNs were randomly placed within the 
ROSAT NEP survey area. In order to reproduce the flux distri- 
bution of the real sample, we followed the method of Mullis 
et al. (2004). In practice since the cumulative AGNs logN-logS 
source counts distribution can be described by a power law, 
S = kS~", with a - 1.3, the differential probability scales 
as 5"*"^^'. Using a transformation method (e.g. Press et al. 
1986, chapter 7) we see that the random flux above a certain 
X-ray flux S um is distributed as 5 = 5/,>„(l -/?)", where p is 
a random number uniformly distributed between and 1 and 
5/,m=2xl0"'^ erg cm"^ s"', i.e. the flux limit of the NEP sur- 
vey. All random AGNs with a flux lower than the flux limit map 
(see Fig. 4 in Henry et al., 2006) at the source position were 
excluded. In order to assign a redshift to these "sources" we 
computed the predicted redshift distribution at the position of 





Fig. 1. Top Panel: The redshift distribution of the AGNs in the 
NEP survey (shaded histogram) and of the randomly gener- 
ated AGNs with the same selection effects (filled histogram). 
Bottom Panel: The redshift distribution of the galaxy clusters 
in the NEP survey (shaded histogram) and of randomly gen- 
erated clusters (filled histogram). The deviations at low z be- 
tween the data and the random sample are mainly caused by 
the NEP superc luster of galaxies (Mullis et al., 2001). 



each accepted source. Once the flux limit at the position where 
the source was randomly placed is known, and denoting with 
<p(L, z) the luminosity function, then the number of sources per 
redshift interval dz is given by 



N(z)dz 



Li ' 



(p(L,z)dV(z)dLdz, 



(3) 



where L,„„, is the minimum luminosity observable at redshift z 
with a local flux limit S;,„, and dV(z) the differential comoving 
volume element. The k-correction does not play any role in the 
calculations since we assumed an average spectral index r=2 
(as in Mullis et al. 2004) for all AGNs. For the luminosity func- 
tion <p(L, z) we took the luminosity-dependent density evolution 
(LDDE) best fit model of Hasinger et al. (2005). A redshift was 
then randomly assigned to each source via Monte Carlo inte- 
gration of the predicted redshift distribution. For galaxy clus- 
ters we applied the same procedure assuming a slope of the 
logN-logS distribution of a- 1 .3, as luminosity function an A-B 
evolving Schechter model (Rosati et al. 2002) with the parame- 
ters obtained by Mullis et al.(2004)]and using a sensitivity map 
specific for NEP clusters ( [Henry et al. 2006). Since for clus- 
ters of galaxies the k-correction is not negligible, we assigned 
to the random clusters an intrinsic spectrum according to a 
MEKAL spectral model with a plasma temperature kT=3 keV 
and a metallicity Z - Q.3Zq. We also applied to L„„„ a "size- 
correction" according to the results of Henry et al. (2006) in 
order to compensate for the missing flux in the X-ray photom- 
etry aperture and the variation of the angular dimensions of the 
object with z. Such a procedure was repeated until we popu- 
lated the survey volume with 37200 random AGNs and 12600 
random clusters (i.e. 200 times more objects than in the real 
data sample). The redshift distribution of the random cluster 
and AGN samples are plotted in Fig. [1] together with the real 
data. To obtain a realistic estimate of the uncertainties of the 
cross correlations we used the bootstrap resampling technique 
described by e.g. Ling, Frenk and Barrow (1986). 

5. Results 

We present the spatial cross-correlation function between clus- 
ters and AGNs in Fig.|2] A positive clustering signal is detected 
in the distance interval s <50 /z^,,' Mpc. In order to test the 
strength of the clustering we performed a canonical power-law 
fit, 



^ca(s) = - 



(4) 



with So and y as free parameters. The fit can be done using 
the coarsely binned data in Fig.|2]and minimizing the;^'^ statis- 
tics. However, this approach is extremely sensitive to the size 
and distribution of the bins. In order to overcome this problem 
we performed a standard maximum likelihood power-law fit to 
unbinned data. The comoving separation s was parsed in very 
small bins so that there are either or 1 data pairs for any given 
interval. In this regime Poisson probabilities are appropriate. 
In order to perform the fit, we need to find the predicted proba- 
bility distribution of the cluster- AGN pairs for each value of y 
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Fig. 2. The Cluster- AGN soft X-ray cross correlation function 
plus one. The error bars are quoted at Icr level. The dashed 
line represents the best fit maximum-likelihood power-law fit 



io=8.7;^P3 /i^y Mpc andy=1.7;^g^. The shaded region illustrates 
the Icr confidence region of the power-law fit in the distance 
range in which it was performed. 



and So- We calculated the number of predicted pairs by replac- 
ing ^ca(s) in Eq. |2] with the model given by Eq. |4]and using 
DcDa(s) (hereinafter A{s)ds) as variable. We can then use the 
separations of all the N cluster- AGN pairs to form a likelihood 
function. This is defined as the product of the probabilities of 
having exactly one pair in the interval ds at each separation 
of the N pairs times the probability of having no pairs in all 
the other differential intervals. This for all 5 in a chosen range 
{sa-Sb), in our case where ^ca{s) can be reasonably represented 




n 



1 10 100 

<s„03a> {h-7o' Mpc) 

Fig. 3. The distribution of the angular radial separations trans- 
lated into redshift space distances at < z >=0.38 



by a power law. Assuming Poisson probabilities we thus obtain 
the likelihood 



X = ]~~[ A{_Si)ds exp 



(5) 



where A{si)ds is the expected number of pairs in the interval 
ds, and the index j runs over all the elements ds which do not 
contain pairs. We then define the usual quantity S = -2 In X 
and drop the terms independent of model parameters (see e.g. 
Schuecker & Boehringer 19981 ICroft et al. 199711 leading to 



;J A(s)ds-2Y^\n[A(sd]. 



(6) 



In order to check in which range of separations we can conduct 
our analysis, we transformed the angular separations between 
clusters and AGNs into an average redshift space separation. 
We assumed that all the sources were at < z >=0.38 (i.e. the 
median redshift of the cluster and AGN sample). The result is 
plotted in Fig. [3] As one can see there are no real pairs with 
separation < 2 h^^ Mpc. For this reason we decided not to con- 
sider points at separation lower than 2.5 /i^q Mpc. These points 
in the cross-correlation function are mainly introduced by the 
parameters RcDa, RaDc in Eq. |2] which could have smaller 
separations than the real data since the random sample includes 
neither the extended emission of the clusters nor the broadening 
due to the PSF of pointlike sources. This prevents us also from 
overestimates of the correlation length introduced by the ampli- 
fication bias (see e.g lVikhhnin & Forman 1995l l. For this rea- 
son and since on scales larger than 50 /z^^q Mpc there is no ev- 
idence of signal, the fit was performed over the distance range 
2.5-50 h^iQ Mpc. In Fig.|4]we show the results of the maximum- 
likelihood power-law fit to ^ca{s) for the ROSAT NEP survey. 
The 1, 2 and 3cr were obtained at A 5 levels of 2.3, 6.2 and 11.8 
from the minimum value of S . The best fit parameters obtained 
are sq-^.I^^q^ h^^^ Mpc and y=1.7^Q^ where the uncertainty 
is at the Icr confidence level. With y fixed to 1.8 (i.e. a typi- 
cal value found in galaxy-galaxy correlation function) we find 
sq ~ 8.5/i^Q Mpc, a similar value was obtained by extending 
the fitting region to 60 hj^ Mpc and restricting it to the 2.5^0 

/l^y Mpc. 

The integral constraint (I Peebles 19801 1. which is a systematic 
shift in correlation functions introduced by the limited vol- 
ume observed, was computed following the prescription of 
[Roche et al. ( 19931] This can be obtained numerically with a 
fit by assuming that the correlation function is represented by 
a power law of fixed index (here we used -y-l.l) on all scales 
sampled by the survey. The underestimate of ^ca{-^) due to the 
integral constraint is found to be of order 2%. 
In order to determine the stability of these results the proce- 
dure was repeated first by separating the field in two subfields 
twice, the North and South, and West and East parts of the sur- 
vey. The fluctuations due to sample variance are found to be 
smaller than the typical amplitude of the uncertainties. A sim- 
ilar result is obtained by recomputing the ^ca('S) for clusters 
withz <0. 18 (i.e. the cluster median redshift) andz > 0.18. The 
dependence of ^ca{s) on the cluster X-ray luminosity (Lx) was 
evaluated dividing the cluster sample into two subsamples with 
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7 a 9 10 11 12 

Sq (h^J Mpc) 

Fig. 4. The Icr, 2cr and 3cr confidence contours in the so,y 
space for the power-law fit to ^ca('?) for two interesting param- 
eters. The cross represents the best fit values. 

Lx < (or >)3.8 X 10'*3 h-^^ erg s '. Though no significant Lv 
dependent behaviour was detected, we cannot yet conclude that 
there is no luminosity or redshift dependent cross-correlation 
length because of the low-statistics of the sub-samples. The 
stability of the result was also checked fitting the data with the 
likelihood estimator used by MuUis et al. (2004) returning no 
significant deviations from our results at 1-2 cr level. 

6. Discussion 

We presented here the first direct evidence of spatial clus- 
tering of soft X-ray selected AGNs around X-ray selected 
clusters of galaxies. Indirect evidence was presented by 
Henry et al. (1991), Cappi et al. (2001), Cappelluti et al. (2005) 
(and references therein). These authors found significant X- 
ray point source overdensities (about a factor 2) around distant 
clusters of galaxies when compared to cluster-free fields. If the 
overdensities were at the cluster redshift they would arise at 
scales smaller than ~7 Mpc. Since the correlation function 

is proportional to (y) , a ^ca = 1 implies an overdensity of a 
factor 2 with respect to a randomly distributed field. We can 
conclude that, since the correlation length found in this work 
reflects the scale of the overdensities known up to now, we ob- 
serve a physical overdensity (of at least a factor 2) of AGNs 
around clusters between 2 and ~ 8 h^^ Mpc from the center of 
the clusters. 

Because of the shallowness of the NEP survey, the AGN sur- 
face density (i.e. <30 deg^ in the central region) does not allow 
detection of such a correlation via overdensity analysis since it 
would be dominated by small number statistics. In fact, from 
our results we expect to detect AGNs overdensities on scales 
<7-8 /i^g Mpc from the center of clusters. At < z >~0.18 (i.e. 
the median z of the cluster sample of the NEP survey) these 
overdensities arise on scales of ~0.6 deg"^ which are easily 
resolved by the NEP survey. However to significantly detect 
these overdensity on single clusters, a conspicuous number of 
sources is necessary to disentangle real overdensities from shot 



noise. This problem could be easily resolved by high angu- 
lar resolution telescopes like Chandra and partially by XMM- 
Newton. In this direction deep and wideQ Chandra and XMM- 
Newton surveys like COSMOS, which will return an AGN sur- 
face density of up to ~ 2700 deg"^, would allow seeing in a 
0.015 deg^ region (i.e. the size of an ACIS-I chip) a population 
of at least 40 AGNs belonging to the cluster environment (i.e. 
assuming an average overdensity of a factor 2, we expect ~40 
sources belonging to the cluster and ~40 to the background). 
As an example, at the limiting flux of the C-COSMOS sur- 
vey, the AGN population of a cluster at z-l would be observed 
with a 0.5-2 keV limiting luminosity L,„,„ ~10'*^ erg s"^ At 
z-\ the size of an ACIS chip (i.e. 8 arcmin) corresponds to a 
linear dimension of ~4 Mpc. Having 40 AGNs in a sphere 
with this radius, corresponds to a space density of AGN with 
Lx >L„„„ of ~0.15 /i^y Mpc"-'. We can however state that, ac- 
cording to the result presented here, clusters of galaxies could 
be detected by AGN overdensities (rather than galaxy overden- 
sities) if the depth of the survey would provide an AGN surface 
density sufficient to overcome the Poisson noise on the AGN 
number In general in order to understand the galaxy evolu- 
tion in dense environments the measure of cross correlation 
between clusters and different kind of galaxies is an impor- 
tant tool. We already know that infrared dusty galaxies avoid 
dense environments therefore showing a large cross-correlation 
length and weak clustering signal in the small separations re- 
gion (Sanchez et al., 2005). We also know that blue galaxies 
avoid low-z rich clusters cores (Butcher & Oemler, 1984). It 
is therefore important to compare the cluster-AGN cross cor- 
relation length to that of cluster and different galaxy type. Mo 
et al. (1993) computed the cross-correlation function of Abell 
clusters and QDOT IRAS galaxies. They found an average cor- 
relation length and a slope in agreement with the results pre- 
sented here. Moreover [Lilje & Efstathiou (1988)| showed that 
the cross-correlation function of Abell clusters with Lick galax- 
ies is positive on scales <29 Mpc with a slope y ~ 2.2 and a 
correlation length of ~ 12.6 /z^q Mpc. These results are also in 
agreement within 1 cr with our findings on AGNs. These first 
comparisons already suggest that AGNs are clustered around 
galaxy clusters just like galaxies. As a final check we compared 
our ^cA to the X-ray cluster-galaxy cross-correlation func- 
tion (hereinafter CGCCF) computed by [Sanchez et al. (2005)| 
They used the X-ray selected clusters of the REFLEX survey 
(Bohringer et al. 2002) and the galaxies from the APM survey 
(IMaddox et al. 1 9901 limited to b;=20.5 mag. They found that 
the CGCCF behaves like a broken power-law with a cut-off dis- 
tance of ~2 hj^^ Mpc with a steeper slope at small distances. We 
can define the following approximate biasing relations: 

^CAis) ^ bcbA^p(s), ^ccis) ^ bcbc^pis). (7) 

Here ^p{s) is the autocorrelation of matter, bc,bA and be are 
the bias factors relative to galaxies, AGNs and clusters, respec- 
tively. By dividing the two equations we can then derive j^(s). 



' The expected fluxes limit for the C-COSMOS and XMM- 
COSMOS survey are 1-2x10"'^ erg cm"^ s"' and -6x10"'* erg cm"^ 
s" , respectively. These surveys will cover 0.9 deg^ and 2.1 deg^ re- 
spectively 
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In order to perform this operation several effects must be taken 
in account. 

- The bias of REFLEX clusters could be slightly differ- 
ent from that of NEP clusters since they are differ- 
ently distributed in redshift and the surveys have dif- 
ferent limiting fluxes and X-ray luminosity distributions. 
[Schuecker et al. (200 1)| showed with the power spectrum 
that even having a large sample of clusters as REFLEX the 
error on the bias determination is still high. Since in the 
NEP survey we expect an even higher uncertainty it is a 
reasonable approximation to consider the ratio consis- 
tent to 1. 

- Sanchez et al. (2005) computed ^cG in real space while 
we work in redshift space. In order to evaluate the effect 
of redshift space distortion on ^cc we used the results of 
[Croft et al. (1999)] Their Figs. 3 and 5 show ^cc computed 
both in redshift and real space. From that work we estimate 
that our relative bias is affected by a ~3Q% overestimate 
below 10 /!7(5 Mpc and of ~10% between 10 h^^^ Mpc and 
20 Mpc. Their work also indicates that ^cc does not de- 
pend on the richness of the clusters used in the calculation 
and that the correction for the scale independent biasing on 
large scales jKaiser 1987 ) can be neglected when compared 
to the size of our uncertainties. 

The ratio t^{s) is plotted in Fig.|5]as a function of the distance 



s {h,„' Mpc) 

Fig. 5. The ratio between the observed ROSAT NEP ^ca(s) and 
the best fit ^cc(r) obtained by [Sanchez et al. (2005)| Errors are 
quoted at the Icr level. The shaded region shows the expected 
level of j^(i)=l if the cross-correlation functions were com- 
pared in the same space. 

7. Summary 

We derived for the first time the soft X-ray spatial cross- 
correlation function between clusters and AGNs using the data 
of the ROSAT NEP survey. A strong clustering signal was de- 
from the center of the cluster The shaded region shows the tected on scale s <50 /i7,J Mpc. The best power-law fit param- 
eters are So=S.7^^'^ /i^q Mpc and ■y-l.l^Q^j. In this work we 
observed that the source density of AGNs is higher near clus- 
ters than in the field. This result confirms earlier findings of 
overdensities of AGNs around clusters reported by many au- 
thors and improves the evidence connecting the overdensities 
to the large scale structure of the Universe. We also derived 
the relative bias between AGNs and galaxies which is consis- 
tent with one on almost all scales investigated here. This result 
still allows, within the errors, a factor 2 fluctuation. New wide 
field surveys (such as XMM-COSMOS) performed with the 
new generation X-ray telescopes will be useful to enlarge the 
statistics, to better understand the physics of AGNs in clusters 
and to extend the analysis to the inner regions of clusters. 



value of our measurement that implies that ^(s)-l when tak 
ing into account the difference between real and redshift space 
measurements discussed in the previous paragraph. The ratio 
is consistent with 1 on almost all scales. We cannot exclude, 
within the errors, much different values of the relative bias. Our 
data suggest an average relative bias consistent with unity but 
allow an upper limit of ~6 (at Icr) at separations s < 50 h^^ 
Mpc. For separations 5 > 10 /z^,,' Mpc no lower limits big- 
ger than zero can be given. On larger scales the error increase 
thus it is difficult to draw any conclusion. At large separations 
the power-law shape of ^cc becomes uncertain, this makes a 
comparison of our data with those of Sanchez et al. (2005) 
less meaningful. We cannot exclude a significant antibiasing 
of AGNs when compared to galaxies, especially at low sepa- 
rations. Though the amplitude of the uncertainties of our data 
still allows a fluctuation in the relative biasing of more than 
a factor 2, we can conclude with a precision of Icr that the 
probability for a galaxy to become an AGN is constant in the 
range of separations sampled in this work and that AGNs can 
be considered as good tracers of the dark matter distribution as 
are galaxies. New deep and wide field surveys such as XMM- 



COSMOS ( |Hasinger et al. 2006 | [Cappelluti et al. 2006) and 



C/ianiira-COSMOS (P.I.: Martin Elvis) with their spectroscop- 
ically identified >2000 AGNs and ~100 clusters will improve 
the precision of this work and will allow us to investigate the 
behaviour of AGNs also in the inner regions of the clusters. 
In this region the observations suggest that the AGNs source 
density is affected by other physical effects such as merg- 



ing ( Ruderman & Ebeling 2005 Branchesi et al. 2006) or ram 
pressure phenomena 
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